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ABSTRACT 

We report high angular-resolution (~1") CO J=3-2 interferometric mapping, 
using the Submillimeter Array (SMA), of IRAS 22036+5306 (122036), a bipolar 
pre-planetary nebula (PPN) with knotty jets discovered in our HST SNAPshot 
survey of young PPNs. In addition, we have obtained supporting lower-resolution 
(~10") CO and 13 CO J=l-0 observations with the Owens Valley Radio Obser- 
vatory (OVRO) interferometer, as well as optical long-slit echelle spectra at the 
Palomar Observatory. The CO J=3-2 observations show the presence of a very 
fast (~220km s^ 1 ), highly collimated, massive (O.O3M ) bipolar outflow with a 
very large scalar momentum (about 10 39 g cm s _1 ), and the characteristic spatio- 
kinematic structure of bow-shocks at the tips of this outflow. The Ha line shows 
an absorption feature blue-shifted from the systemic velocity by ~100km s -1 , 
which most likely arises in neutral interface material between the fast outflow 
and the dense walls of the bipolar lobes at low latitudes. Since the expansion 
age of the outflow as determined from our data is only about 25 years, much 
smaller than the time required by radiation pressure to accelerate the observed 
bipolar outflow to its current speed (^ 10 5 yr), the fast outflow in 122036, as in 
most PPNs, cannot be driven by radiation pressure. The total molecular mass 
(0.065 Mq) is much less than that derived from a previous detailed model of the 
near to far-infrared SED of this object (showing the presence of a large, cool 
dust shell of total mass 4.7M ), most likely due to the interferometric observa- 
tions resolving out the CO flux from the dust shell due to its size. We find an 
unresolved source of submillimeter (and millimeter-wave) continuum emission in 
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122036, implying a very substantial mass (0.02-0.04 M & ) of large (radius^l mm), 
cold (^50 K) dust grains associated with I22036's toroidal waist. We also find 
that the 13 C/ 12 C ratio in 122036 is very high (0.16), close to the maximum value 
achieved in equilibrium CNO-nucleosynthesis (0.33). The combination of the high 
circumstellar mass (i.e., in the extended dust shell and the torus) and the high 
13 C/ 12 C ratio in 122036 provides strong support for this object having evolved 
from a massive (^4M ) progenitor in which hot-bottom-burning has occurred. 

Subject headings: circumstellar matter - planetary nebulae: individual (IRAS 
22036+5306) - reflection nebulae - stars: AGB and post-AGB - stars: mass loss 
- stars: winds, outflows 



1. INTRODUCTION 

Pre-Planetary nebulae (PPNs) - objects in transition between the AGB and planetary 
nebula (PN) evolutionary phases - hold the key to one of the most vexing and long-standing 
problems in our understanding of these very late stages of evolution for low and intermediate 
mass stars. Current observational evidence shows that during the PPN phase, the dense, 
isotropic mass-loss which generally marks the late AGB phase changes to high-velocity, bipo- 
lar or multipolar mass-loss. The emergence of fast collimated outflows or jets has been hy- 
pothesized as the primary mechanism for this dramatic change in the geometry and dynamics 
of the mass-loss (Sahai & Trauger 1998). However, the physical mechanism for producing the 
fast outflows remains unknown. (Sub)millimeter-wave interferometric observations offer one 
of the best probes of the dynamics and energetics of the shock acceleration process, which 
transfers a substantial amount of directed momentum to large parts of the dense AGB wind. 

IRAS 22036+5306 (122036) is a bipolar PPN with knotty jets discovered in our HST 
SNAPshot survey of young PPNs, with an estimated distance and luminosity of 2kpc and 
23OOL , respectively (Sahai et al. 2003 [SZSCM03]). Our HST imaging suggested that 
122036 is a key object for understanding how jets can sculpt bipolar lobes in a progenitor 
AGB star wind, since its jets and their working surfaces appear to be simultaneously ob- 
served. In this paper (the second one by Sahai and collaborators on this important PPN) 
we report CO J=3-2 observations of the PPN with the Submillimeter Array (SMA), which 
resolve a massive high-velocity bipolar outflow aligned with the optical nebulosity and clearly 
pinpoint the location and nature of the bow-shocks at the tips of this outflow. We also report 
supporting lower-resolution (~10") CO and 13 CO J=l observations with the Owens Valley 
Radio Observatory (OVRO) interferometer, as well as optical long-slit echelle spectra at the 
Palomar Observatory. 
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The rest of the paper is organised as follows: in §2 we describe our observational 
techniques and data-reduction, in § 3 we describe our observational results, in § 4 we discuss 
the kinematics, age and mass of the bipolar outflow, in § 5 we derive the total molecular 
mass and compare it to the mass of dust in the circumstellar environment, in §6 we derive 
and discuss the implications of the 13 C/ 12 C isotope ratio in 122036, in § 7 we discuss the 
formation of the nebula, and in § 8 we summarize our main conclusions. 

2. OBSERVATIONS 

The CO J=3-2 observations were obtained with the SMA, using two tracks, on July 
20 (UT) and Aug. 29 (UT), 2004, with the array in the compact and extended configura- 
tions, respectively. The J2000 coordinates used for 122036 were RA = 22 fe 05 m 30.28 s , Dec 
= +53°21'33.0". The receivers were tuned to center the 345.796 GHz rest frequency of the 
CO (3-2) line in the upper sideband (USB), after correcting for the source's -45 kms -1 LSR 
velocity. On the night of July 20th, the weather was excellent, with 0.7 mm of precipitable 
water vapor. The antennas' DSB system temperatures ranged from ~ 200K near the zenith 
to ~ 350K at an elevation of 30°. Unfortunately one of the eight antennas had extremely 
high phase noise on that night, and its data was unusable. On the night of Aug. 29, all 
eight antennas were working well, but the precipitable water vapor varied between 1.0 to 1.6 
mm, causing the system temperatures to range from ~ 250K to ~ 500K. For both tracks 
the correlator was configured to have the highest available uniform resolution, 406.25 kHz 
channel (0.35 kms -1 ), over a total bandwidth of 2 GHz. Observations of Saturn and Uranus 
were obtained for bandpass calibration, and Titan and Uranus were used for flux calibration. 
The quasars BL Lac and 2013+370 were observed every 1/2 hour for use as complex gain 
calibrators. A total of 12 hours of on-source data was taken, split nearly evenly between the 
two tracks. 

The data were calibrated using the MIR package, and mapping was done with AIPS. 
Visibility weights were derived from receiver Tsys measurements. After calibration, the 
data was exported to AIPS for mapping. Two spectral lines, CO (3-2) and an unexpected 
second line, were detected in the USB, but no lines were seen in the LSB. A continuum level 
was fitted to the regions of the USB data which were free of line emission, and subtracted 
from the visibility data. Two data cubes were then generated. When producing the first 
cube, we included data from both tracks with the aim of generating maps of the emission 
integrated over large velocity intervals with the smallest synthesized beam. For the second, 
we discarded the more noisy data from the extended track, with the aim of generating high- 
quality emission maps with much narrower velocity bins. The final data cubes had 4872 
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channels after the removal of overlapping portions of the spectra lying in the rolloff regions 
of the hybrid correlator's analog filters. In the first (second) data cube, the clean beam 
had a FWHM= l'/12 x l'.'OO (1"45 x 1"29) at a position angle (PA) of 88°.6 (-84°. 7). A 
continuum map was made using the lower sideband (LSB) data, and uniform weighting of 
Tsys-weighted visibilities. The central 250 MHz of the LSB IF coverage was rejected, in 
order to guard against contamination from the strong CO emission in the USB. A 180 kA 
taper was applied to the visibility data when maps were produced. Uniform weighting was 
used. 

Interferometric observations of the CO and 13 CO J=l-0 line emission were obtained 
at OVRO on 2002 September 21 as part of a snapshot survey of pre-planetary nebulae. 
General details on the observational setup are given in Sanchez Contreras & Sahai (2004); 
data reduction and calibration steps are similar to those described in Sahai et al. (2005). 
The clean beam had a FWHM of 9"8 x 5"9 (PA 87°) for CO and 10'.'2 x 6"7 (PA 82°) for 
13 CO. Scale conversion factors were 1.6 and 1.5 K/(Jy beam" 1 ) for CO and 13 CO. 

Since 1 22036 displays a prominent Ha emission line, and the presence of a fast neutral 
outflow can produce a self-absorption signature in this line (e.g., Sanchez Contreras & Sahai 
2001), we also obtained a high-resolution (AA=0.4A) spectrum of 122036 in order to resolve 
the velocity structure of the line, using the Palomar 60-inch spectrograph in echelle mode. 
In this mode, the spectrograph provided a total wavelength coverage of -3600-9000A, with 
a resolution of about 19000. The echelle grating had 52.65 lines/mm and a blaze angle of 
63°26'. Cross-dispersion was provided by two fused quartz prisms with a 60° apex angle. The 
slit was 7'.' 4 long, and had a width of l'/4. The spectra were recorded on a thinned, lumogen 
coated TI 800 x 800 chip. A Thorium-Argon lamp was used for wavelength calibration, and 
the standard star HR7596 was used for flux calibration. The data were obtained during the 
nights of 2002, August 4 and 5, and reduced and calibrated using standard IRAF techniques. 



3. RESULTS 

3.1. Spectral Profiles 

The (spatially) averaged CO J=3-2 line profile shows a strong central core and weaker 
extended wings (Fig. la). The wings cover a total velocity extent (FWZI) of ~450 km s _1 and 
are most likely associated with a fast post-AGB outflow. The large width (FWHM^70 km s _1 ) 
of the central core of the CO J=3-2 emission indicates that much of the core emission is also 
associated with the fast outflow. 

Like the J=3-2 line profile, the CO and 13 CO J=l-0 line profiles obtained also show 
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wide wings (^100 km s -1 ) but do not cover as wide a velocity range (Fig. lb). The observed 
extent of the CO J=l-0 line wings is limited by the spectrometer bandwidth (^200 km s -1 ), 
whereas that of the 13 CO line wings is sensitivity limited. The J=l-0 line profile has a 
triangular shape with a total velocity extent (FWZI) of greater than 100 km s^ 1 , with no 
obvious central core. It is likely that this J=l-0 emission is associated with the material 
seen in the central core of the J=3-2 line. The 13 CO J=l-0 line shows a narrow central 
core, with a width of about 20 km s -1 , and weaker wings extending over a similar velocity 
range as the CO J=l-0 line. Since the typical expansion velocities of AGB circumstellar 
envelopes (CSEs) are 10-15 km s -1 , the narrow central core seen in the 13 CO J=l-0 line may 
be associated with the remnant progenitor AGB envelope in 122036, or slowly expanding 
gas in the dense waist. 

In addition to the CO J=3-2 line, we detected a line at a rest frequency of 346.53 GHz, 
most likely due to the SO N=8-7, J=9-8 transition with a frequency of 346.5284810 GHz 
(JPL catalog). The observed SO line is much narrower than the CO line, with a full width 
of 21km s -1 , and an integrated flux of 20Jykm s _1 (Fig. lc). The SO emission appears to 
be slightly elongated and extended over ~ 2", with a PA of about —30° (Fig. 2, top-left 
corner panel). Within the relatively low signal/noise, the red and blue wings of the line are 
co-located, indicating that the SO emission does not arise in a bipolar outflow. Given its 
similar velocity-width compared to the narrow core of the 13 CO 1-0 line, the SO emission is 
most likely associated with the same material as the former. 

The Ha line shows a profile with a central peak, very wide wings, and an absorption dip 
in the blue wing of the line, centered at Vlsr of about —145km s -1 (Fig. Id). This peculiar 
Ha line shape is seen in many PPNs and young, compact PNs (van de Steene, Wood, & 
van Hoof 2000, Sanchez Contreras & Sahai 2001, Sahai & Sanchez Contreras 2004, Sahai 
et al. 2005). The blue and red wings can be seen extending to velocities of roughly —400 
and 500km s -1 , respectively. The Ha line profile (e.g. IRAS 19024+0044, Hen3-1475) is 
most simply interpreted in terms of the following model - the (spectrally) broad Ha emission 
arises in a compact, spatially-unresolved central source 1 ; this emission and stellar continuum 
is scattered by dust in the extended, expanding nebula. The blue-shifted absorption arises 
due to a largely atomic outflow within the foreground lobe, which absorbs the scattered 
photons. The scattering is expected to produce an overall red-shift of the Ha line peak 
relative to the systemic radial velocity, as observed - the peak of the Ha line is located 
at Vlsr ~ km s -1 , red-shifted from the systemic velocity, Vlsr = —45 km s -1 . Such a 
model has been successfully used to quantitatively explain the spatio-kinematic distribution 



1 the nature of this source is unknown, although possible mechanisms include emission from a very high 
velocity outflow, Raman scattering and/or Keplerian rotation in a dense disk 
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of similar Ha emission profiles observed with STIS/HST in the PPN Hen 3-1475 (Sanchez 
Contreras & Sahai 2001). 

3.2. Velocity-Channel Maps 

We show the spatial distribution of the CO J=3-2 emission as a function of velocity 
in Fig. 2. The centroid of the emission is seen to shift steadily from west to east in the 
panels as the average velocity decreases from 217 to -217 km s -1 , indicating a strong velocity 
gradient along the major axis of the nebula. A plot of the central ±15km s" 1 range with 
narrower velocity bins of width 2.47 km s _1 (Fig. 3) also shows roughly a East- West velocity 
gradient. We cannot separate a kinematically distinct component in the J=3-2 emission in 
the velocity range covered by the central core of the 13 CO J=l-0 line. The overall velocity 
gradient can be seen more easily in a plot of the blue and red wing emission of the CO 
J=3-2 line, overlaid on the optical image of 122036: there is a large separation between 
the red and blue-shifted emission, implying the presence of a bipolar outflow (Fig. 4) along 
the long axis of 1 22036. In order to register the HST and SMA images, we first rotated 
the latter to the orientation of the former. For the final fine translational registration, we 
assumed that peak "C" SZSCM03 coincides with the peak of the integrated CO J=3-2 
intensity (shown in Fig. 5) since the absolute astrometry of the HST images, relative to that 
of the SMA images, cannot be determined to a sufncently high accuracy (i.e., better than 
fewxOVl). The outflow is tilted such that the W-lobe (E-lobe) is on the near (far) side - this 
tilt is consistent with that inferred from the difference in brightness between the two optical 
lobes and the morphology of the waist (SZSCM03). We have not shown the corresponding 
integrated intensity or velocity channels maps for the OVRO observations, because these do 
not resolve the CO J=l-0 emission due to the large beamwidth. 

The map of the emission integrated over all velocities Fig. 5 shows that the the emission 
lies in an elongated region, oriented along the major axis of the optical nebula. A central, 
unresolved, continuum source (at A=0.88mm) with a flux of 0.29±0.04 Jy is also seen; the 
error is dominated by calibration uncertainties. The peak of the continuum source appears 
to be offset from the line intensity peak by ~0'/23. An unresolved continuum source is also 
detected in the OVRO data, with a flux of 8.4± 1.7 mJy 
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3.3. Position- Velocity Structure 

A position-velocity plot along the bipolar axis (Fig. 6) shows abrupt increases in radial 
velocity at offsets of about ±1" from the center. Such structure is typical of what is expected 
when a high-velocity bipolar collimated outflow or jet (hereafter simply referred to as the 
jet) interacts with dense ambient material, producing bow-shocks at the head of the jet 
(e.g., see models and Fig.2, right panel, in Lee & Sahai 2003: hereafter LS03). The bow- 
shocks implied by the CO velocity structure appear to be located in the general vicinity of, 
and may be identified with, the bright regions SI, S2 in the W-lobe, and their presumed 
optically obscured counterparts in the E-lobe. The P-V structure has a high degree of 
symmetry around the center, indicating that the bipolar outflow and surrounding medium 
are correspondingly symmetric. 

In addition to the high-velocity bow-shock structures, the P-V diagram shows a weaker 
and significantly more spatially-extended component, which extends out to the tips of the 
optical bipolar nebula. This component is composed of relatively low-velocity emission, but 
the limits of its velocity structure at the high end cannot be separated from the bow-shock 
emission. 

There are two possible explanations for the simultaneous presence of the low- and high- 
velocity emission structures, with (1) their different spatial extents seen in the CO J=3-2 P-V 
plot and (2) correspondence to the morphological features seen in the optical images. One 
possibility is that the extended lobes in 122036 are geometrically distinct from the outflow 
which produces the high-velocity emission and the bright optical knots S1,S2. The latter is 
probably a collimated outflow whose axis is tilted significantly towards the observer, whereas 
the extended lobes are more in the plane of the sky (assuming that the 3-D outflow velocity 
of the material in these lobes is not too different from that in the high- velocity outflow). 
So in 3-D, 122036 may be a quadrupolar nebula, with the near-side lobe representing the 
S1,S2 outflow seen in projection in front of the extended lobe. A second possibility is that 
the extended lobes have a similar inclination as the material in the fast bipolar outflow, and 
the emission with low (projected) velocity in the P-V plot comes from an outflow which is 
intrinsically slower than the former. 

A position- velocity plot along the minor axis of the 1 22036 shows no significant gradient 
(Fig. 7). 
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4. The Post-AGB Bipolar Outflow 

4.1. Kinematics and Expansion Age 

We have evidence of fast outflows in 122036 from three independent probes. Two of 
these are the CO and Ha data in this paper, and a third is the VLA mapping of OH maser 
emission reported in SZSCM3 and Zijlstra et al. (2001). The OH maser features cover a total 
velocity spread of about 76 km s _1 , and are roughly aligned along the nebular axis, with a 
spatial extent comparable to the distance between the CO-identified bow-shock locations. 
The P-V signature of the OH is thus similar to that of CO except at the bow-shock locations, 
where there is an absence of OH emission extending to the large outflow velocities seen in 
CO. Either OH is destroyed, or the appropriate conditions for exciting observable maser 
emission are not met, in and near the bow-shock region. 

Both the OH and CO emission presumably come from molecular gas in the dense walls 
of an elongated lobe presumably resulting from the interaction of a jet with an ambient 
circumstellar envelope as shown by LS03. In contrast to OH and CO emission, the Ha 
absorption probes atomic gas. Since the Ha absorption feature is not spatially resolved in 
our data, it is difficult to definitively establish its relationship to the fast molecular outflow. 
In the jet-envelope interaction model, the atomic gas may be located in the interior of 
the lobes, representing the pristine, unshocked material of the underlying jet, and/or the 
interface between the latter and the lobe walls. We would thus expect the observed velocity 
of this atomic material to be generally larger than that of the dense material everywhere in 
the lobe. But the blue-shift of the Ha absorption feature from the systemic velocity implies 
a projected outflow speed of about 100 km s -1 , which is smaller than that seen in the CO 
at the bow-shocks. This discrepancy can be explained in the context of the LS03 models, if 
the Ha absorption feature is produced in interface material (but not "pristine", unshocked 
jet material) at low latitudes in the outflow lobe, since the velocity of the interface material 
near the lobe walls is found to decrease with latitude (where the outflow axis represents the 
polar axis) whereas that of the unshocked jet is constant. 

The LS03 modelling shows that, in the momentum-driven case for a jet-envelope in- 
teraction, the velocity field of much of the dense material in the resulting elongated lobe 
is predominantly radial; the age of the outflow is well approximated by dividing the radial 
distance of the bow-shock at the lobe-head by the radial velocity of the material there. If we 
take 220 km s _1 and 1", respectively, as the projected radial velocity and radial distance of 
the material at the bow shock, based on Fig. 6, and if we assume an intermediate inclination 
angle (i = 30°) of the outflow axis to the sky-plane, then the fast outflow age is only 25 (tan 
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V0.58)yr. 2 

4.2. High- Velocity Outflow: Mass and Scalar Momentum 

We now estimate the mass, Mu P , and scalar momentum, P sc in the high-velocity bipolar 
outflow from the integrated J=3-2 line profile shown in Fig. la. In this analysis, we have 
used the J=3-2 emission flux over a velocity range extending from 424 (-505) km s _1 to 15 
(-15) km s -1 from the systemic velocity for the red (blue) component of the outflow. We 
have thus excluded the emission in a velocity range ±15 km s" 1 centered on the systemic 
velocity, which, as suggested by the 13 CO J=l-0 and the SO line profiles, might come from 
a kinematic component distinct from the bipolar outflow. We find values of the velocity- 
integrated J=3-2 emission flux over the blue and red components of the outflow to be 168.3 
and 192.5 Jykm s" 1 , repectively. 

4-2.1. Excitation Temperatures and Optical Depths 

We need to determine the excitation temperature and evaluate any optical depth effects 
in order to compute the mass of emitting material. The ratio of the integrated fluxes of the 
the CO J=l-0 and 3-2 lines, R13, over the high- velocity range covered by both the OVRO 
and SMA observations (i.e. an interval of 215 km s _1 centered at the systemic velocity, but 
excluding the central ±15 km s _1 ) is 0.081, lower than the value in the central ±15 km s _1 , 
which is 0.12. For optically thin emission in both J=l-0 and 2-1, R13 = 0.081 implies an 
excitation temperature, T ex =lA.7K and an emitting mass of 0.022 M Q . However, if either or 
both of lines are not optically thin, then the inferred temperature and mass will be different. 

In order to accurately estimate optical depths as a function of velocity, a detailed spatio- 
kinematic model is needed, however, we do not attempt that here since the emission is 
not sufficiently well-resolved in our data, and so the detailed geometrical structure of the 
emission is poorly known. Nevertheless, rough estimates of the optical depths can be made by 
considering the sizes, total fluxes and total velocity dispersion of the low- and high-velocity 
components. We first deconvolve the J=3-2 line profile into a central gaussian with a width 
(FWHM) of 30 km s" 1 which represents the low- velocity component, and two pairs of high- 



2 The much longer estimate of the age (250 yr for an inclination of 30°) given by SZSCM03 was due to 
(i) the use of a factor 5 smaller outflow speed based on the OH data, and (ii) it applied to knot K\ which is 
located a radius about twice as large as that of the bow shocks 
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velocity gaussian components (HVlr,HVlb and HV2r,HV2b) whose sum represents most of 
the emission from the high-velocity component. Assuming that the centers of the near (far) 
pair, HVlr and HVlb (HV2r and HV2b), are symmetrically offset from the systemic velocity 
by ±V1 (±V2), and their width is Wl (W2), we find the values of VI, V2, Wl, W2, as well 
as the peak intensities of the gaussians which provide the best fit to the J=3-2 line. We 
find Wl~92km s _1 and W2~220km s -1 . The total flux in the low- velocity component is 
100 Jy-km s _1 , whereas HV1 and HV2 have about 191 and 170 Jy-km s -1 each. 

The angular extent of the high-velocity component is estimated to be about 1"2, i.e., 
3.6 xlO 16 cm from the J=3-2 blue and red wing emission maps (Fig. 4). In this estimation, we 
have used geometric means of the major and minor axis of both the emission regions and the 
beam to obtain deconvolved sizes. Although the detailed spatio-kinematic structure of the 
emitting regions is not known, if we assume that the large velocity gradient approximation 
is valid, we can obtain rough estimates of the optical depths. We have simulated a large 
velocity gradient by assuming a spherically expanding cloud with constant expansion velocity 
equal to half the FWHM of the model components determined above. Such a cloud has 
significant velocity gradients along most lines-of-sight (except those which pass very close to 
the center). For each of the four high- velocity gaussian components above, we have computed 
the optical depths (and thus the masses) and excitation temperatures needed to fit the R13 
ratio and produce half the CO J=3-2 flux of each component (since the fractional flux within 
the FWHM points of the 2-dimensional gaussians representing each of the HV1 and HV2 
components is 0.5). These masses are then scaled up by a factor of 2 in order to obtain the 
mass corresponding to the total flux in each of the HV components. 

We have assumed a fractional /co=CO/H 2 abundance of 2 x 10~ 4 (SZSCM03). For each 
of the red and blue parts of the HV1 (HV2) component, we find T ex =18 (15.5) K, line-of-sight 
J=3-2 optical depth at a radius half the radius of the model cloud, r 32 of 1.1 (0.35), and a 
mass of 4.4 x 10~ 3 (3.1 x 10~ 3 ) M Q . Adding together the HV1 and HV2 masses from the red 
and blue parts, and doubling the resulting value to account for the mass contributing to the 
flux outside the half-power points of the gaussian fits to these ccomponents, we get a total 
emitting mass in the high-velocity outflow of M bip =0.03 M . 

In order to calculate the scalar momentum, we have used the formulation as described 
in Bujarrabal et al. (2001), assuming optically-thin emission, and a single temperature for 
the high- velocity outflow, which we take to be the emission-weighted average (17 K) of the 
HV1 and HV2 components above. We find P sc — 1.1 x 10 39 g cm s _1 for an intermediate 
inclination angle (i = 30°) of the nebular axis to the sky-plane. This outflow cannot be driven 
by radiation pressure: the dynamical time-scale of the lobes, for % = 30°, is only ~25yr, 
much smaller than the time required by radiation pressure to accelerate the observed bipolar 
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outflow to its current speed, t ra d ^ P S c/{L/c) ~ 10 5 yr, given I22036's luminosity of 23OOL . 
The value of t ra( i is a lower limit because only a fraction of the total luminosity, namely that 
directed into the solid angle representing the opening angle of the bipolar outflow, may be 
utilised for driving the latter. 

5. Molecular Mass versus Dust Shell Mass 

The size of the low-velocity component is estimated to be about l'/3 (4xl0 16 ), from the 
J=3-2 emission map in Fig. 3. The flux and angular size of this component are comparable 
to that in each of the red and blue components in the HV1 and HV2 features, but it is spread 
over a much smaller velocity range. Hence it is expected to be significantly more optically 
thick than the HV1 and HV2 components, making the mass determination more uncertain. 
Fitting this component as above, we can find reasonable fits with T ex ~ 26 — 33 K, r 32 about 
7.5-10, and a mass of about 0.042 — 0.035 M . The corresponding optical depth of the J=l-0 
line, Tio is ~2.6-2.2 in these models, hence any significant reduction in the total mass of the 
shell reduces r w below unity in large parts of the emitting region, resulting in a value of 
R13 which is too low compared to the observed value of 0.12. For example, if we reduce Tio 
by a factor 2, we get R13=0.08, inconsistent with our data. Thus the minimum mass of the 
low-velocity component is about 0.035 M . 

Our data support the relatively high values of T ex for the low-velocity component derived 
from our modelling. The relatively high surface brightness observed at the systemic velocity 
(2.2 Jy/beam, Fig. 3) implies a brightness temperature of about 12 K. Accounting for a 
beam-dilution factor of 2 (since the source size is about the beam size), the actual source 
brightness temperature is T;,=24K, and the corresponding excitation temperature required 
for optically thick emission is 32 K. 

Thus the total molecular mass is M tot =0.065 M . This is much lower than the value 
(~4.7M ) estimated from the dust mass (~O.O23M ) inferred from a model of I22036's near 
to far-infrared fluxes (as observed by ISO, MSX and IRAS) (SZSCM03), assuming a gas-to- 
dust ratio, 5=200. In this model, cool (67-35 K) dust resides in a large spherical shell (with 
inner and outer radii of 1.4 x 10 17 and 5.3 x 10 17 cm) of radial optical depth Ay ~ 1.5, and 
hot (1000-350 K) dust resides in a very optically-thick (A v ~ 120), compact ((3-4.6) x 10 14 
cm) central torus. The mass of the cool shell is 4.7M , and the (very poorly constrained) 
mass of the compact central torus is 6 x 1O _4 M , assuming a 100/im dust emissivity of 34 
cm 2 g _1 . An additional warm (145-82 K) shell (with inner and outer radii of 1.3 x 10 16 and 
5.7 x 10 16 cm) with relatively little mass ~O.OO5M is added to the first two components to 
obtain a good fit to the spectral energy distribution (SED) in the ~20-40 /im range. 
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The size of the CO J=3-2 emission source along its long axis is ~ 5", which is signifi- 
cantly smaller than the central cavity (of size ~ 9") in the model cool dust shell. Thus the 
observed CO J=3-2 emission does not probe the material in the cool shell, which explains 
the discrepancy between the mass estimate based on the CO data (~0.032M Q ) and that 
derived from the dust emission (4.7M Q ). The size of the dust shell is large (inner and outer 
radii of 4'/ 6 and 17'/6) compared to the typical size of structures which would not be largely 
resolved out by the interferometer, which is ~ 9" given our minimum UV spacing of 10 
klamba, and good UV coverage for spacings greater or equal to about 23 klambda. Hence it 
is likely that a significant fraction of the CO J=3-2 flux from this shell is resolved out by our 
interferometric observations. Furthermore, astrophysical mechanisms for the lack of observ- 
able CO J=3-2 emission from the dust shell are also plausible, e.g. (a) photodissociation 
of CO due to the interstellar ultraviolet radiation, and (b) insufficient excitation due to low 
gas kinetic temperatures. We think it is unlikely that photodissociation is important, since 
the photodissociation radius (estimated using Table 3 from Mamon et al. 1988) given the 
high mass- loss rate inferred for the shell (5.5 x 1O _4 M yr _1 : SZSCM03), and an assumed 
expansion velocity of 15 km s _1 , is a few x 10 18 cm, which is significantly larger than the 
outer radius of the dust shell. As in the case of other PPNs (e.g., IRAS 19024+0044: Sahai 
et al. 2005, Roberts 22: Sahai et al. 1999), it is plausible that the gas kinetic temperatures 
in the extended dust shell of 1 22036 are quite low because of the combined effect of cooling 
due to adiabatic expansion and the lack of dust frictional heating. 

5.1. The Central Submillimeter Continuum Source: Large Grains in I22036's 

Dusty Torus 

The 0.88 mm and 2.6 mm fluxes reported here are significantly in excess of the SED 
resulting from the dust model of SZSCM03. This can be seen in in Fig. 8, which shows 
a modified version of the SZSCM03 model (discussed below) - the SMA and OVRO data 
points lie significantly above the dashed red curve representing emission from the cool shell, 
which is the dominant component at long wavelengths in the SZSCM03 model. Since the 
observed 0.88 mm continuum source is unresolved, its size is less than 0'.'9 x 0'.'7 (the FWHM 
of the SMA clean beam). Thus the 0.88 mm flux in the SZSCM03 model, within this small 
aperture, would be even lower than shown by the dashed red curve in Fig. 8. Since the 
spectral index of the continuum from 0.88 mm to 2.6 mm, a (defined as F u ~ u a ) is 3.3, 
too steep to be a result of thermal brehmstrahhlung, the submillimeter and millimeter-wave 
fluxes in 122036 are most likely the result of thermal emission from a population of quite 
large, cold grains. We believe that the grains must be large because the opacity spectral 



-13- 



index j3 = din k/ din v in the submillimeter-millimeter wavelength region from our data 3 is 
1.3, which is significantly smaller than the value expected for standard small grains (f3 Ri2). 
Such small values of (3 are consistent with those obtained in theoretical models with grain 
sizes (in a power-law size distribution) extending to the millimeter range (Draine 2006). 
Given the compact nature of the 0.88 mm continuum source, it is quite likely that these 
large grains are physically located within the dusty waist seen in the HST images. Note that 
the bulk of the grains producing the SED of 1 22036 in the ~few-200 /im range are probably 
small (~0.1 /im) as has been found for other PPNs from the very large fractional polarisations 
typically observed towards these objects at optical and near-infrared wavelengths (e.g. Sahai 
et al. 1998a, 1999, Oppenheimer et al. 2005). 

In order to provide rough constraints on the physical parameters of the large-grain dust 
source, we have used DUSTY (as in SZSCM03) to compute the emission from a dust shell 
with large grains, varying its major inputs - the temperature at the shell inner radius (the 
ratio of the outer to inner radius is arbitrarily assumed to be 2), the optical depth at 100 /im, 
r(100), the grain temperature at the inner radius T^(in), and the grain size, a - to fit the 
observed SMA and OVRO continuum fluxes. We find that dust grains with a ^1 mm and 
Td{in) ^ 50 K are needed for producing the observed fluxes at 0.88 and 2.6 mm. The shell 
has an optical depth r(100) = 0.1. Models with smaller grains (e.g. a <^ 500 jim) which fit 
the 0.88 and 2.6 mm fluxes produce too much flux in the ~100-200 /jm range. In order for the 
large-grain shell to be smaller than or equal to the observed size of the 0.88 mm continuum 
source, i.e. < 2.5 x 10 16 cm, the source luminosity corresponding to the required irradiating 
intensity must be ^35OL , i.e., ~15% of the total source luminosity. This result requires 
that the large grains be confined to a region where they can be substantially shielded from 
the full radiation of the central star, supporting our conjecture that they are located within 
the toroidal waist seen in the HST images of 1 22036, since this region is partially shielded 
from the full stellar radiation field by the compact, optically-thick central dust component 
inferred by SZSCM03. From our model, we estimate the mass in dust of this large grain 
component to be ~O.O4M , taking the dust mass extinction coefficient, k, at A << a (say 
at 100 /im) to be 2.5 cm 2 g _1 , derived from k = a(47ra 3 prf/3) _1 , with a = na 2 , a=l mm and 
p d = 3 gem" 3 . 

Our dust model is based on commonly-used values for the wavelength-dependent com- 
plex refractive index of silicate dust particles, however uncertainties in the dust mass ab- 
sorption coefficient may allow for a reduction of the mass estimate by as much as a factor 
of 2. Thus it is clear that there is a very substantial mass of such grains in the dusty waist 



3 derived assuming optically-thin emission in the Rayleigh-Jeans limit, for which a = 2 + (3 
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of 122036. If we assume a gas/dust ratio of 200 in the waist, as in the cool dust shell, the 
total mass in the toroid becomes 4-8M , a disturbingly large value, especially when com- 
bined with the 4.7 Mq of material in the cool dust shell. We therefore conclude that the gas 
content in the waist must be significantly depleted, and the gas-to-dust ratio consequently 
much lower than the typical value, maybe by as much as an order of magnitude. 

6. The 13 C to 12 C isotope ratio 

We derive the 13 C to 12 C abundance ratio, /13/12, in 122036 from the ratio of the J=l-0 
13 CO and 12 CO fluxes. The difference in the shapes of the 12 CO and 13 CO lines shown in 
Fig. 1 is consistent with a uniform value of /13/12, but larger optical depths in the central 
component compared to the high- velocity component. As indicated by the prominent core 
emission seen in the 13 CO but not the 12 CO line profile, the ratio of the 13 CO to 12 CO flux in 
the core (0.37) is much higher than that in the high- velocity wings (0.2). This difference in 
the flux ratio between the low-velocity and high-velocity outflows is due to the difference in 
their optical depths: t w is about 2 in the central component but 0.4 in the HV1 component 
(see §5). We have therefore computed /13/12 from the flux ratio in the more optically-thin 
part of the lines - the 13 CO/ 12 CO flux ratio of 0.2 gives /13/12 = 0.16. We find that the same 
value of /13/12 reproduces the observed 13 CO to 12 CO flux ratio in the core if Tio=2.2 and 
T ex = 30 K, supporting our modelling of the 12 CO emission from the low- velocity outflow, 

Thus the 13 C to 12 C abundance ratio in 122036 is close to the maximum value achieved 
in equilibrium CNO-nucleosynthesis (0.33), and provides a strong test of models of stellar 
nucleosynthesis. For comparison, other evolved objects which show high values of /13/12 (i.e., 
in excess of 0.1) are (i) the circumstellar envelopes of some oxygen-rich AGB or post-AGB 
objects such as OH231.8+4.2 (Sanchez Contreras et al. 2000), Frosty Leo (Sahai et al. 2000), 
U Equ (Geballe et al. 2005), (ii) the rare class of J-type carbon stars (e.g., Y CVn or T Lyr: 
Lambert et al 1986), and (iii) selected planetary nebulae (e.g., Ml-16: Sahai et al. 1994). 
In general, there is increasing observational evidence that s large fraction of evolved stars 
have 12 C/ 13 C abundance ratios lower than about 15, significantly lower than those predicted 
by standard models. For example, from millimeter-wave CO studies, Palla et al. (2000) and 
Balser, McMullen & Wilson (2002) have estimated the /13/12 ratio in modest PN samples 
(16 & 9 objects, respectively) and find it lies typically in the range ~0.04-0.1. 

There are two processes which can enhance the /13/12 ratio above the values pre- 
dicted by the standard model. The first, cold-bottom-processing (CBP), is a deep mixing 
mechanism believed to operate during the RGB phase of stars with mass less than about 
2M Q ( Wasserburg, Boothroyd & Sackmann 1995), and the degree of enhancement decreases 
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with increasing stellar mass. In a 1M Q star, CBP can produce a value of /13/12 as high as 
0.2 during the RGB phase, but this ratio is lowered to 0.025-0.05 during the AGB phase due 
to the addition of fresh 12 C. Since the progenitor central star of 122036 is significantly more 
massive than 2M & considering the high mass of its dust shell, CBP cannot account for its 
high /13/12 ratio. 

The second mechanism, hot-bottom-burning (HBB), is believed to operate in AGB stars 
with initial masses above about 3.5-4M Q . Such stars develop deep convective envelopes with 
very high temperatures at the base causing the operation of the CN cycle (i.e., 12 C — > 13 C 
— > 14 N). In this case, /13/12 will be near the equilibrium ratio of 0.33, and the surface /13/12 
is expected to asymptotically approach this limit. HBB is thus the likely mechanism for 
producing the high /13/12 ratio in 122036. Rubin et al. (2004) reach a similar conclusion for 
the PN2440, where they find /13/12 ~0.2 from analysis of CIII lines in IUE data. 



7. DISCUSSION 

7.1. The shaping of bipolar Planetary Nebulae 

The dynamical evolution and shaping of PPNs is believed to result from the shock 
interaction between a fast, collimated post- AGB wind and the slowly expanding, dense 
wind ejected during the previous AGB phase. The bipolar post-AGB wind impinges on, 
and thus accelerates and heats, two localized regions of the (presumably spherical) AGB 
wind. A significant amount of momentum is believed to be transferred to the rest of the 
AGB shell by means of bow-like shocks. Although morphological identifications of bow- 
shocks have been made in PPNs, their characteristic spatio-kinematic structure has seldom 
been identified. Thus, the bow-shock structures inferred kinematically from our CO J=3-2 
imaging of 1 22036 provide support for, and are an important input for quantitative modelling 
of the physics of, the interaction of collimated fast winds with AGB envelopes as a production 
mechanism for bipolar PPNs (see, e.g., LS03). 

We have shown that the bipolar molecular outflow in 1 22036 cannot be driven by radia- 
tion pressure. Bujarrabal et al. (2001) infer the presence of fast molecular outflows in a large 
sample of PPNs, from the extended wings seen in their CO line profiles, and conclude that 
these outflows cannot be driven by radiation pressure. Only for a small fraction of objects 
in their sample (e.g., CRL2688 or the Egg Nebula: Cox et al. 2001, IRAS 09371+1212 or 
the Frosty Leo Nebula: Castro-Carrizo et al. 2005, CRL618: Sanchez Contreras et al. 2004, 
OH231.8+4.2: Alcolea et al. 2001) has interferometric mapping been carried out, directly 
revealing the structure of the fast outflow. In all such well-mapped PPNs, the fast outflow or 
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each of its components (if it is comprised of multiple components) appears to be collimated 
and directed along well-defined axes. Our results show that 1 22036 fits well in this picture of 
the fast molecular outflows in PPNs. Given the increasing direct evidence for the presence of 
jet-like outflows in PPNs (e.g., W43A: Imai et al. 2002; IRAS 16342-3814: Sahai et al. 2005, 
IRAS 19134+2131: Imai et al. 2004), we believe that the case for the fast molecular outflows 
in 122036 and other PPNs being driven by collimated jets is quite strong. The presence 
of such jets, operating in an episodic fashion and changing their directionality, has been 
proposed as a general mechanism for explaining the diverse morphologies and widespread 
presence of point-symmetry in planetary nebulae (Sahai and Trauger 1998, Sahai 2004). The 
launching mechanism for such jets remains unknown, but given the similarity between the 
fast outflow speeds in PPNs and jets in low-mass pre-main sequence stars with accretion 
disks, it is plausible that the mechanisms which drive jets in PPNs and YSOs are very sim- 
ilar. In this scenario, the accretion disk in a PPN would be produced around a relatively 
close companion low-mass star, which accretes material from the dense AGB wind of the 
primary (Morris 1987, Mastrodemos & Morris 1998). 

Our CO J=3-2 observations and the dust model of SZSCM03 indicate an apparent 
absence of dense cool material in 1 22036 along non-polar directions, in the volume between 
the molecular outflow and the inner radius of the cool dust shell (hereafter, the "in-between" 
region). In the SZSCM03 model, the inferred inner radius of this dust shell cannot be 
significantly smaller than the derived value because I22036's SED (which peaks near 60 fj,m) 
dictates that the emitting grains have low temperatures, and therefore be located relatively 
far from the illuminating source. Even if we explicitly account for the reduction in the 
illuminating flux due to the dusty torus in SZSCM03's model, which covers ~0.4 of the 47r 
solid angle around the central star, the lower limit to the fraction of the central star's total 
luminosity which irradiates the cool dust shell is 0.6. Hence the inner radius of the shell 
(which scales as the square-root of the illuminating flux) is at most 23% smaller than derived 
by SZSCM03. Although one may explain the absence of CO emission using insufficient 
excitation as an argument, the lack of thermal dust emission from the in-between region is 
strong evidence for the lack of dense material there. 

A faint halo has been found in 122036 at 0.8 /^m in new deep images taken with the 
HST/ACS camera (Sahai et al. 2006, in preparation) which indicates the presence of a rather 
tenuous surrounding remnant AGB envelope. But we do not know if the fast molecular 
outflow in 122036 was produced by a collimated jet interacting with this tenuous envelope. 
If the ambient material is of sufficiently low density, then the swept-up shell can be too faint 
to be observable, and one can then see the jet material directly, as appears to be the case for 
the young PN, Hen 2-90 (Sahai & Nyman 2000, Lee & Sahai 2004). Only detailed numerical 
simulations of the jet- AGB envelope interaction in 1 22036 can tell us if the material in its 
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fast molecular outflow (a) represents the swept-up shell due to the interaction of the jet 
with the tenuous AGB envelope seen in the near-IR halo, (b) is intrinsic, i.e., belongs to the 
underlying jet itself, or (c) represents material associated with the dense torus that has been 
accelerated by the jet. In model c, one can imagine that the jet punched its way out of a 
flattened dusty cocoon around the central star as envisaged for the case of the well-studied 
PPN, CRL2688, Sahai et al. (1998b) - cocoon matter entrained by the jet produces the 
observed bipolar outflow in CO emission, and the remnant cocoon represents the torus. 

While considering the above issues, it is useful to note that the nebular properties of 
122036 are very similar to that of the well-studied object, OH231.8+4.2, although their cen- 
tral stars have different spectral types (~F5 versus M9, respectively). As in 122036, OH231.8 
harbors a massive molecular bipolar outflow oriented along the long axis of the optical bipo- 
lar nebula, and a weak round near-IR halo is seen around its bipolar outflow, signifying the 
presence of a tenuous spherical envelope surrounding this outflow (i.e., characterised by a 
mass-loss rate of ~ 1O~ 6 M yr _1 : Alcolea et al. 2001). Furthermore, 122036 also has a very 
high 13 C to 12 C ratio, similar to that in OH231.8. Thus it appears that both 122036 and 
OH231.8 are very similar objects and are formed by a similar mechanism; the main difference 
between the two is that OH231.8+4.2 is a more extreme version - in it, the mass of cool 
material in the bipolar outflow is estimated to be significantly larger (Alcolea et al. 2001). 
Both objects show high- velocity OH maser emission. Since OH231.8 was not observed by 
ISO, we do not know its far-IR SED as well as that of 122036, however the existing IRAS 
data show that the SEDs of both objects are very similar. For example, in OH231.8 (as 
in 122036), (a) the broad-band IRAS fluxes from 12 to 100 /xm show a peak at 60 /im, and 
(b) the LRS data shows a broad silicate absorption feature at 10 /im and a steeply rising 
spectrum towards longer wavelengths. OH231.8 also shows a compact continuum component 
at 1.3mm (Sanchez Contreras et al. 1998), which the authors interpret as emission from 
grains larger than about 5 /im. 

122036 most likely had a fairly massive progenitor, easily ^4M . The combination of 
the high circumstellar mass (i.e., in the extended dust shell and the torus) and the high 
13 C/ 12 C ratio in this object provides strong support for the operation of the hot-bottom- 
burning mechanism. 

7.2. Dusty Waists and Large Grains 

Compared to the question of forming bipolar and multipolar morphologies in PPNs, 
the formation of the dense, dusty waists in these objects remains a much bigger puzzle. 
The presence of large grains with a very substantial mass associated with 1 22036's toroidal 
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waist, is a surprising, but not unprecedented result. In the well-studied PPN, AFGL2688, 
Jura et al. (2000) find evidence for a substantial mass (~O.O1M ) of large-sized grains. In 
the Red Rectangle (RR: also the nearest PPN to us at a distance of 350-700 pc, bipolar in 
morphology and harboring a disk in rotation around a central binary), Jura et al. (1997) 
derive a minimum mass of 5 x 1O~ 4 M in large (a ^0.02 cm), cool (50 K) grains based on 
RR's 1.3 mm continuum flux. Men'shchikov et al. (2002) fit the full SED and near-infrared 
images of the RR using a 2-D dust radiative transfer model, and infer a dust mass of about 
O.O1M in grains with a ~1 mm in a compact, central torus. In the quadrupolar PPN 
IRAS 19475+3119, Sarkar & Sahai (2006) find, from a comparison of a detail model of its 
full SED (from optical wavelengths to 200/zm) to the observed 0.85 mm flux, that there is 
a substantial 0.85 mm excess flux which also requires the presence of large, cool grains. We 
speculate that the formation of dusty waists and the presence of large grains in PPNs are 
due to intimately-linked physical processes, thus a resolution of the issue of waist formation 
may well lead lead to a resolution of how the large grains are made or vice versa. A possible 
mechanism is the destruction of volatile comet ary- debris disks within several hundred to a 
thousand AU (Stern et al. 1990) by an intermediate-mass star during its luminous post- 
main-sequence evolution, providing a source of large solid particles which could form a dusty 
torus. However, the total mass of large grains which we have found in 1 22036 is far too large 
compared to the mass estimates for our Kuiper Belt (~0.1M e , e.g., Luu & Jewitt 2002) or 
the inner parts of the Oort Cloud (~40M®, e.g., Dones et al. 2004), and argues against such 
a process, unless of course the Kuiper Belt/Oort Cloud analogs of an intermediate mass star 
like that in 1 22036 are significantly more massive than those in our Solar system. 



8. Conclusions 

Using the SMA, we have made high angular resolution (~ 1") maps of the CO J=3-2 
emission in IRAS 22036+5306, a bipolar pre-planetary nebula with knotty jets, discovered 
previously in one of our HST imaging surveys of candidate PPNs. These data have been 
supplemented by lower-resolution CO and 13 CO J=l-0 observations with OVRO, as well as 
optical long-slit echelle spectroscopy. We find 

(i) The CO J=3-2 spatially integrated line profile covers a total velocity extent (FWZI) of 
~500km s -1 . The CO and 13 CO J=l-0 also show large velocity widths, although with a 
smaller spread due to the lower sensitivity of these data. 

(ii) The CO J=3-2 emission comes from a massive high-velocity bipolar outflow aligned 
with the optical nebulosity. The characteristic position-velocity signature of a bow-shock 
resulting from a high- velocity bipolar collimated outflow/jet interacting with dense ambient 
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material can be seen in the CO J=3-2 data in both the red and blue-shifted components of 
the outflow. The P-V structure has a high degree of symmetry around the center, indicat- 
ing that the bipolar outflow and surrounding medium are correspondingly symmetric. The 
bow-shocks do not lie at the tips of the optical lobes. 

(hi) The Ha line shows a profile with a central peak, very wide wings, and an absorption 
dip in the blue wing of the line, blue-shifted from the systemic velocity by about 100 km s -1 . 
The absorption feature indicates the presence of a fast, neutral, outflow inside the lobes. 

(iv) A central, unresolved, continuum source with a flux of 0.29±0.04Jy (8.4±1.7mJy) at 
A=0.88mm (2.6 mm) is seen in the SMA (OVRO) data. 

(v) The mass in the high- velocity outflow, i.e., at velocities offset more than ± 15 km s _1 from 
the systemic velocity, is 0.03 M , and the scalar momentum is l.lxlO 39 g cm s _1 . The large 
value of the scalar momentum, together with the small age of the high-velocity outflow, 
implies that it cannot be driven by radiation pressure. 

(vi) The total molecular mass derived from integrating the flux of the CO J=3-2 line over its 
full velocity range, is ~O.O65M , much lower than the value (~4.7M ) estimated previously 
from a dust-shell model of 1 22036's near to far-infrared fluxes, assuming a gas-to-dust ratio, 
5=200. This is probably due to the interferometric J=3-2 observations resolving out the 
flux from the large shell (inner and outer radii of 1.4 x 10 17 and 5.3 x 10 17 cm) in which the 
cool dust resides. Furthermore, the gas in this shell may not be sufficiently excited. 

(vii) The small size of the continuum source indicates that it is associated with the dusty 
toroidal waist of 122036. The submillimeter and millimeter-wave continuum fluxes imply 
the presence of a very substantial mass (~O.O2-O.O4M ) of large (of radius ^ 1 mm) cold 
grains in this waist. 

(viii) The 13 C/ 12 C ratio in 122036 is 0.16, close to the maximum value value achieved in 
equilibrium CNO-nucleosynthesis (0.33). The enhancement is most likely due to the oper- 
ation of hot-bottom-burning, which operates in stars with initial masses larger than about 
4M . 

We thank Nick Sterling for pointing out the importance of hot-bottom-burning in en- 
hancing the 13 C/ 12 C ratio, and Noam Soker for reading an earlier version of this paper. RS 
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Table 1: IRAS22036+5306 Model Components 



Component 


Parameter 


Value 


Comments 


Cool, 


Rjra 


1.4xl0 17 cm 


SZSCM03 


Dense Shell 


Rout 


5.3xl0 17 cm 






T,(in) 


67K 






Td(out) 


35 K 






density law 








Mass 


4.7M 






M 


5.5xlO~ 4 M yr" 1 






grain size 


MRN 




Warm Shell 




1.3xl0 16 cm 


SZSCM03 




Rout 


1.9xl0 16 cm 






T d (in) 


145 K 






T d (out) 


82 K 






Mass 


5.5xlO~ 3 M 






grain size 


MRN 




Inner Torus/Disk 


Rjra 


3xl0 14 cm 


SZSCM03 




Rout 


4.5xl0 14 cm 






T d (in) 


1000 K 






T d (out) 


365 K 






Mass 


6xlO- 4 M 






grain size 


MRN 




Toroid 


Rjra 


9xl0 15 cm 


this paper 




Rout 


1.8xl0 16 cm 






T d (in) 


50 K 






T d (out) 


35 K 






Mass(dust) 


O.O2-O.O4M 






grain size 


a ^1 mm 




Fast Mol. Outflow 


Speed 


-220 km s- 1 


this paper 




Age 


25 yr 






Mass 


O.O3OM 






Momentum 


9.2 xlO 38 g cm s" 1 






trad 


^ 10 5 yr 




Slow Mol. Outflow 


Speed 


—15 km s _1 


this paper 




Mass 


O.O35M 




Fast Atomic Outflow 


Speed 


-100 km s- 1 


this paper 
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Fig. 1. — The spectra of the (a) CO J=3-2 line integrated over the entire emission region, (b) 
CO and 13 COl-0 peak intensity in the unresolved emission source for each line, (c) the SO 
9(8)-8(7) line integrated within a 2"x2" area centered on the SO peak (5-channel Hanning 
smoothing has been applied), and (d) the Ha line (flux in erg s _1 cm _2 A _1 normalised to a 
peak value of unity), in IRAS 22036+5306 
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Fig. 2.— Maps of the CO J=3-2 and SO 9(8)-8(7) emission in 122036 obtained with the 
SMA. For SO, the emission is quite weak and only the (velocity) integrated intensity is shown 
(top left corner). The contours on the SO map are -30, 30, 60 and 90% of the peak flux 
of 0.54 Jy/beam. For the CO J=3-2 line, the remaining panels show emission maps as a 
function of velocity. Each panel covers a velocity bin of width 9.86 km s _1 , centered at the 
velocity offset (as measured from the systemic velocity Vlsr = —45 km s" 1 ) shown in the 
top right corner. The contour levels are -5, 5, 10, .. 95% of the peak flux in the central 
channel which is 2.2 Jy/beam. The center of each panel (marked by a cross) is located at 
RA = 22 A 05 m 30.28 s , Dec = +53°21'33.0". The la noise is about 38 mJy/beam 
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Fig. 3. — As in Fig. 2, but with the total velocity coverage extending over the central 
±15km s _1 range, with each panel showing emission in a velocity bin of width 2.47km s -1 . 
Contour levels are -10,10,20,. ..90% of the peak flux in the central channel which is 2.2 
Jy/beam 
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Fig. 4. — Map of the blue and red wing CO J=3-2 line emission from IRAS 22036+5306, 
overlaid on a 0.6 /xm HST image. The velocity ranges used to define the blue and red wings 
in this plot are offset by -245 to -61km s _1 , and 61km s _1 to 245 km s _1 from the systemic 
velocity of -45 km s _1 . The peak surface brightness in the red (blue) lobe is 4.10 (3.37) Jy- 
km s _1 /beam. Countours are 0.82, 1.23, 1.64... 3.28 Jy-km s _1 /beam in the blue lobe, (i.e. 
with a step of 0.41 Jy-km s _1 /beam), and 0.82, 1.23, 1.64, ... 3.69, 3.95 Jy-km s _1 /beam in 
the red lobe. The size and orientation of the clean beam (1'.'12 x l'.'OO, PA 88°. 6) are shown 
in the lower-right corner 
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Fig. 5. — Maps of the velocity-integrated CO J=3-2 line emission [magenta contours; levels 
are 15 to 95% of the peak emission (13.3 Jy-km s _1 /beam) in steps of 10%] and 0.88 mm 
continuum emission [green contours; levels are 15 to 95% of the peak emission (0.17 Jy/beam) 
in steps of 10%] from IRAS 22036+5306, obtained with the SMA, overlaid on a 0.6 /an HST 
image (processed to enhance sharp features, as in SCZCM3). The la noise in the continuum 
map is 4.5 m Jy/beam. The sizes and orientations of the clean beams for the line maps 
(1'.'12 x l'.'OO, PA 88°. 6) and continuum maps (0"94 x 0"72, PA 83°. 9) are shown in the 
lower-right corner 
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Fig. 6. — Position- velocity plot of the CO J=3-2 emission along the major axis of the nebula. 
Contour levels are 5 to 95% of the peak in steps of 5%. A color HST image of the nebula 
(with the 0.6 /im image in green and the 0.8 /zm image in red, as in SZSCM03) is shown, 
with its major axis aligned along the spatial axis in the P-V plot, for comparison 
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Fig. 7. — Position- velocity plot of the CO J=3-2 emission along the minor axis of the nebula. 
Contour levels are 10,20. ..90% of the peak in steps of 10% 



31 




Log (Wavelength [microns]) 

Fig. 8. — Observations [black curve: ISO spectra, blue symbols: photometric data - SMA 
0.88 mm and OVRO 2.6 mm: circles w/crosses, MSX: squares, IRAS: triangles, ISO/PHOT: 
cross] and a model spectrum (magenta curve) of 122036. Individual components (dashed 
curves) of the model are also shown: cool (red) & warm (cyan) shells, hot inner disk (blue) 
and a component with large, cold grains (green), probably associated with I22036's dusty 
waist. The flattened region at the long-wavelength end of the ISO spectrum is most likely 
due to contamination from interstellar cirrus emission. 



